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1.  Introduction 


Since  the  discovery  of  the  detonation  [1,2],  the  question  of  physical  mechanisms  that  create  this  self-supporting, 
supersonic,  shock-driven  reaction  wave  has  been  at  the  forefront  of  the  combustion  theory.  The  development  of  a 
detonation  is  a  significant  threat  to  chemical  storage  and  processing  facilities,  mining  operations,  etc.  [3],  while  con¬ 
trolled  detonation  initiation  in  the  next  generation  of  propulsion  systems  promises  to  revolutionize  transportation  [4], 
On  astrophysical  scales,  detonation  formation  is  presently  believed  to  be  the  most  important,  yet  least  understood, 
aspect  of  the  explosion  mechanism  [5,  6]  powering  type  la  supernovae  (SNIa),  the  thermonuclear  incineration  of  a 
degenerate  white  dwarf  star.  The  use  of  SNIa  as  cosmological  standard  distance  indicators  has  led  to  the  discovery  of 
the  accelerating  expansion  of  the  Universe  [7,  8],  suggesting  the  existence  of  dark  energy.  Future  high-precision  cos¬ 
mology  studies  of  dark  energy  will  require  accurate  calibration  of  SNIa  as  standard  candles,  which  will  be  impossible 
without  understanding  the  process  of  detonation  formation. 

Already  first  systematic  studies  of  detonation  showed  [2]  that  it  can  arise  from  a  slow,  highly  subsonic  deflagra¬ 
tion  ignited  in  an  initially  unpressurized  system.  Significant  progress  has  since  been  made  experimentally  [9,  10] 
and  numerically  [11,  12,  13,  14,  15]  in  elucidating  the  physics  of  the  deflagration-to-detonation  transition  (DDT)  in 
confined  systems,  in  particular  in  closed  channels.  These  studies  showed  that  the  confining  effect  of  channel  walls  on 
the  hot,  expanding  burning  products  and  the  interaction  of  the  resulting  flow  with  walls  and  obstacles  are  instrumental 
in  providing  significant  flame  acceleration  and  pressure  increase,  thus  creating  the  conditions  necessary  for  detonation 
ignition.  This  raises  the  question:  is  DDT  possible  in  unconfined  media  without  assistance  of  walls  or  obstacles,  e.g., 
in  unconfined  clouds  of  fuel  vapor  or  in  the  interior  of  a  white  dwarf  during  a  SNIa  explosion? 

It  was  originally  suggested  by  Zel’dovich  et  al.  [16]  that  a  detonation  can  form  in  a  region  (hot  spot)  with  a  suitable 
gradient  of  reactivity.  The  resulting  spontaneous  reaction  wave  propagating  through  that  gradient  creates  a  pressure 
wave  that  can  eventually  develop  into  a  shock  and  a  detonation  [17,  18].  In  confined  systems,  multidimensional  direct 
numerical  simulations  (DNS)  have  shown  that  hot  spots  can  form  through  repeated  shock-flame  interactions  and  fuel 
compression  by  shocks  [11], 

It  remains  unclear,  however,  if  and  how  hot  spots  would  form  in  unconfined,  unpressurized  media.  The  most 
likely  mechanism  is  flame  interaction  with  intense  turbulence.  In  particular,  it  was  suggested  [18,  19]  that  disruption 
of  the  internal  flame  structure  by  high-speed  turbulence,  producing  a  distributed  mode  of  burning,  can  create  hot 
spots  capable  of  initiating  a  detonation.  It  is  unknown,  however,  whether  this  can  indeed  occur,  as  there  are  no 
realistic  ab  initio  experimental  or  numerical  demonstrations  of  this  process.  Here  we  show  that  high-speed  turbulence- 
flame  interaction  can  indeed  lead  to  DDT,  but  through  a  different  process  that  does  not  rely  on  the  propagation  of  a 
spontaneous  reaction  wave  and,  thus,  does  not  require  the  formation  of  hot  spots. 

2.  Physical  model  and  numerical  method 

The  DNS  calculations  presented  here  solve  compressible  reactive-flow  equations  with  thermal  conduction,  molec¬ 
ular  species  diffusion,  and  energy  release  [20,  21],  They  use  an  ideal  gas  equation  of  state  and  first-order  Arrhenius 
kinetics  that  describes  chemical  reactions  converting  fuel  into  product.  A  simplified  reaction-diffusion  model  rep¬ 
resents  stoichiometric  H2-air  and  CH4-air  mixtures  under  Le  -  1  conditions  with  model  parameters  calibrated  to 
reproduce  both  laminar  flame  and  detonation  properties  [12,  13].  The  reaction  model  for  CH4-air  was  validated 
against  experimental  data  on  DDT  in  obstructed  channels  [13]. 

Simulations  were  performed  with  the  code  Athena-RFX  [20,  21,  22],  It  uses  a  fully  unsplit  corner  transport 
upwind  scheme  with  PPM  spatial  reconstruction  and  the  HLLC  Riemann  solver  [23]  (see  [22]  for  the  detailed  tests 
of  the  hydrodynamic  solver  and  [20,  21]  for  the  discussion  of  the  reactive-flow  extensions).  Turbulence  driving  is 
implemented  via  a  spectral  method  [24,  20] . 

3.  Numerical  simulations 

We  consider  the  interaction  of  a  premixed  H2-air  flame  with  high-speed,  steadily  driven  turbulence.  The  sim¬ 
ulation  setup  is  similar  to  our  previous  studies  [20,  21],  which  analyzed  the  quasi-steady  turbulent  flame  evolution 
(simulation  6  in  Fig.  1).  Here  we  consider  a  larger  system  and  a  higher  turbulent  intensity  (simulation  10  in  Fig.  1). 
The  computational  domain  is  a  uniform  256  x  256  x  4096  Cartesian  mesh  with  width  L  =  0.5 1 8  cm,  giving  the  resolu¬ 
tion  Ax  =  ()/(,/ 1 6,  where  d/o  ~  0.032  cm  is  the  laminar  flame  thermal  width  in  cold  fuel.  Kinetic  energy  is  injected  at 
the  scale  L  to  produce  a  homogeneous,  isotropic  turbulence  with  characteristic  velocity  U  =  1.9  x  104  cm/s  «  6 3. S'  Lq 
at  the  scale  L,  where  S l.o  =  3.02  X  102  cm/s  is  the  laminar  flame  speed  in  cold  fuel.  The  corresponding  large-scale 
eddy  turnover  time  is  Ted  —  L/U  =  27.3  /rs,  while  integral  velocity  (7/  =  1.2  X  104  cm/s  ~  405 l.o  and  integral  scale 
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Figure  1:  (Color  online)  Combustion  regime  diagram  [25]  showing  performed  DNS  calculations.  Symbol  color  and  shape  indicate  the  reactive 
mixture  considered  and  the  observed  mode  of  burning.  The  full  flame  width  Ip, o  ~  [20]. 


I  —  0.12  cm.  Resulting  turbulence  away  from  the  flame  has  equilibrium  Kolmogorov  energy  spectmm  <x  k  5/  !  in  the 
inertial  range  extending  to  scales  <  b/o  [20], 

Atf  =  0,  fuel  is  at  the  temperature  To  =  293  K  and  pressure  Pq  =  1.01xl06erg/cm3.  The  planar  flame  is  initialized 
normal  to  the  longest  dimension  of  the  domain  (z-axis)  with  zero-order  extrapolation  z-boundary  conditions  and 
periodic  boundaries  along  other  directions.  After  ~2reci ,  the  turbulent  flame  becomes  fully  developed,  and  it  reaches 
a  quasi-steady  state  (QSS)  which  lasts  until  t  ~  6.5 Telj.  Figure  2  shows  the  turbulent  flame  speed,  S  j,  based  on  the 
fuel  consumption  rate  [20],  Turbulent  flame  properties  during  this  period  are  consistent  with  the  previous  analysis 
of  such  QSS.  In  particular,  the  flame  remains  in  the  thin  reaction  zone  regime  [20]  with  the  reaction  zone  structure 
virtually  unaffected  by  intense  turbulence.  S  j  is  primarily  controlled  by  the  increase  of  the  flame  surface  area  with 
the  additional  occasional  increase  <30  -  40%  due  to  periodic  flame  collisions  and  the  formation  of  cusps  [21]. 

In  contrast  with  the  behavior  at  a  smaller  JJ\  and  in  a  smaller  system  (case  6),  here  the  QSS  lasts  a  relatively  short 
time  (Fig.  2),  and  after  t  ~  6.5 red  the  flame  evolution  changes  qualitatively:  S  j  increases  rapidly,  becoming  supersonic 
by  7. 1 8rr(/  and  exceeding  the  Chapman-Jouguet  detonation  velocity,  Dqj,  at  1.5tcij.  DDT  occurs  shortly  thereafter  at 
7.53 Tej,  and  S  j  reaches  its  maximum  at  7.58^.  At  7.63rfy/  a  fully  developed  overdriven  planar  detonation  emerges, 
and  it  quickly  relaxes  to  Dcj. 

The  system  evolution  during  this  process  is  shown  in  Fig.  3.  At  6.39rc</  a  slight  overpressure  has  formed  inside  the 
flame  brush,  and  the  energy  generation  rate  per  unit  volume,  E,  is  still  close  to  its  value  in  the  planar  laminar  flame. 
As  the  pressure  grows  and  the  turbulent  flame  accelerates,  fuel  inside  the  flame  brush  is  compressed  and  heated.  This 
increases  the  local  flame  speed,  S causing  E  to  rise.  Note,  at  later  times  E  exceeds  its  laminar  value  by  ~2  orders 
of  magnitude.  Such  accelerated  burning  leads  to  further  fuel  compression  and  larger  S l-  The  resulting  feedback 
loop  drives  a  catastrophic  runaway  process  that  produces  a  large  pressure  build-up  and  creates  strong  shocks.  Before 
such  shocks  are  able  to  form  a  single  global  shock,  their  collision  creates  a  high-pressure  triple  point  that  ignites  a 
detonation  (details  of  this  last  stage  will  be  presented  in  a  separate  paper). 

To  determine  the  regime  of  burning  during  the  runaway,  we  recorded  the  average  temperature,  T /,  and  pressure, 
Pf,  of  pure  fuel  (Y  >  0.95)  inside  the  flame  brush.  Up  until  the  moment  of  DDT,  T /  remains  <  700  K,  and  the 
corresponding  induction  times  are  significantly  larger  than  all  dynamical  timescales.  Furthermore,  at  all  times,  the 
average  internal  flame  structure  reconstructed  using  method  described  in  [20]  is  close  to  that  of  a  laminar  flame  in  fuel 
with  the  corresponding  T /  and  Pf.  Thus,  during  the  runaway,  burning  is  controlled  by  flame  propagation  and  not  by 
autoignition,  which  excludes  the  possibility  of  formation  of  global  spontaneous  reaction  waves. 

4.  Mechanism  of  the  spontaneous  runaway 

Consider  an  unconfined  fluid  volume  V  with  the  total  internal  energy  s.  To  increase  the  pressure  inside  V  (as  in 
Fig.  3),  an  energetic  process  must  generate  energy  comparable  to  s  on  a  characteristic  sound-crossing  time  of  this 
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Figure  2:  (Color  online)  Turbulent  flame  speed,  St,  normalized  by  the  sound  speed  in  cold  fuel,  cSto.  Legend  gives  simulation  numbers  (Fig.  1). 
Shaded  gray  regions  show  the  range  of  critical  values  of  St  according  to  eq.  (1)  based  on  the  sound  speed  in  fuel,  csj,  and  product,  cs^p,  for  fuel 
temperatures  in  the  range  320  -  430  K.  Blue  dots  on  the  curve  for  simulation  10  indicate  times  of  individual  profiles  in  Fig.  3.  Time  is  normalized 
by  the  corresponding  value  of  reci  in  each  simulation. 


volume,  i.e.,  s  ~  s/ts.  If  this  volume  represents  a  flame  with  width  5  and  cross-sectional  area  L2,  i.e.,  V  =  5L 2,  then 
the  burning  speed  of  such  flame  is  defined  as  S  —  m/pfL2  [20],  where  m  =  s/q  is  the  total  fuel  consumption  rate 
and  pf  is  the  fuel  density.  Then  condition  s  ~  e/ts  can  be  rewritten  as  S  ~  csE/qpf ,  where  ts  =  6/cs,  cs  is  sound 
speed,  and  E  =  e/V  is  the  internal  energy  per  unit  volume.  The  flame  here  may  be  laminar,  turbulent,  or  distributed, 
provided  it  has  the  required  burning  speed. 

To  show  the  physical  meaning  of  this  condition  on  S ,  assume  ideal  gas  equation  of  state,  E  =  P/(y  -  1).  At 
the  start  of  the  runaway,  pressure  is  nearly  constant  across  the  flame.  Then  product  density  is  pp  =  pfTf/Tp  = 
p/T f/(T f  +  q/Cp)  =  P/(P/pf  +  q(y-  l)/y),  where  Tp  is  the  product  temperature  and  Cp  is  the  specific  heat  at  constant 
pressure.  For  energetic  reactive  mixtures,  the  denominator  P/pj  +  q(y  -  1  )/y  can  be  approximated  as  q(y  -  1).  For 
instance,  in  our  case,  q  =  43.28 RTq/M  »  Po/po-  and  at  the  onset  of  the  runaway  P  «  (1.5  -  2)Pq  and  pf  «  po,  giving 
the  accuracy  of  this  approximation  ~  6%  -  11%.  Thus,  pp  ~  P/q(y  -  1),  and  we  finally  get 


cs  cs  P  cs 

— E=—-f - -*-=SCj 

qpf  pf  q(y  -  1)  a 


(1) 


where  a  -  Pf/pp  is  the  fluid  expansion  factor. 

In  the  reference  frame  of  a  steady  flame,  ppUp  =  p/U /  =  p/S ,  where  Uf  and  Up  are  the  velocities,  respectively, 
of  fuel  and  product.  Thus,  eq.  (1)  is  equivalent  to  the  statement  that  Up  =  cs.  If  cs  is  taken  to  be  the  sound  speed  in 
product  then  the  flame  satisfying  eq.  (1)  is  a  Chapman- Jouguet  (CJ)  deflagration  [26], 

The  speed  of  a  CJ  deflagration,  S  cj,  is  a  theoretical  maximum  for  the  flame  burning  speed.  The  discussion  above 
shows  that  such  a  flame  generates  enough  energy  on  its  sound-crossing  time  to  raise  its  internal  pressure  and,  thus, 
disrupt  its  steady-state  structure.  Real  laminar  flames,  both  chemical  [26]  and  thermonuclear  [18],  do  not  have  burning 
speeds  that  approach  S  cj ■  Turbulent  flames,  however,  can  develop  such  high  values  of  S  j. 

Unlike  a  laminar  flame,  in  which  local  sound  speed  smoothly  increases  from  its  value  in  fuel,  cs,/,  to  that  in 
product,  cs,p,  a  turbulent  flame  effectively  consists  of  two  fluids  with  either  csj  or  cs,p.  Thus,  to  verify  that  eq.  (1) 
is  indeed  the  criterion  for  the  onset  of  runaway,  we  show  in  Fig.  2  S  cj  based  on  both  csj  and  cs,p-  Fuel  heating  by 
turbulence  causes  csj  and  cs,p  to  increase  and  a  to  decrease.  Shaded  gray  areas  show,  for  both  sound  speeds,  the  range 
of  critical  values  of  S  j  corresponding  to  fuel  temperatures  320  -  430  K.  In  particular,  in  simulation  10,  T f  w  320  K 
at  2Ted  (lower  bound  of  the  shaded  regions)  and  it  increases  to  ~430  K  by  6.5  r,Y/  (upper  bound). 

Figure  2  shows  that  upon  first  reaching  the  QSS,  St  is  close  to,  but  still  below,  csj/a ,  which  prevents  the  onset 
of  the  runaway.  During  the  time  (2  -  6.5 )red,  however,  fuel  heating  by  turbulence  causes  S  l  to  increase  by  a  factor  of 
~  2,  thus  accelerating  S  t  above  the  critical  value  csj/a  and  allowing  the  runaway  to  begin.  Figure  3(c)  shows  that  at 
this  point  the  product  velocity  indeed  becomes  ~csj.  Furthermore,  growth  rate  of  S  j  increases  significantly  once  S  t 
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Figure  3:  (Color  online)  The  x-y- averaged  profiles  of  (a)  pressure,  P,  (b)  fuel  mass  fraction,  Y,  (c)  z- velocity,  uz,  and  (d)  energy  generation  rate  per 
unit  volume,  E,  in  simulation  10.  Time  since  ignition  for  each  profile  is  shown  in  panel  (b)  (these  times  are  indicated  with  blue  dots  in  Fig.  2).  The 
E  is  normalized  by  its  value  in  a  planar  laminar  flame  propagating  in  cold  fuel,  E^o  =  qS  L,opo/dL,o,  where  q  is  chemical  energy  release  and  po  is 
cold  fuel  density. 


becomes  >cs^p/a,  i.e.,  when  U p  becomes  supersonic  relative  to  both  sound  speeds.  Note  also  that  the  transition  from 
a  QSS  to  a  detonation  occurs,  effectively,  on  a  sound-crossing  time  of  the  turbulent  flame  ts  =  St/c^o  ~  27/rs  »  Terj, 
where  £>7-  »  1  cm  is  the  flame  brush  width  (Fig.  3b)  and  c^o  ~  3.7  x  104  cm/s. 

For  comparison.  Fig.  2  also  shows  S  j  for  turbulent  FF-air  flames  for  other  values  of  (7/  and  l  (Fig.  1).  In  particular, 
in  simulations  5-7,  S  j  remains  well  below  csj/a  and  the  flame  evolves  in  the  QSS  similar  to  that  described  in  [20,  21]. 
This  QSS  was  observed  over  significantly  longer  periods  of  time  than  shown  in  Fig.  2,  e.g.,  16 Teci  in  case  6.  Cases  1-4 
had  similar  behavior  and,  thus,  are  not  shown.  We  did  observe  the  runaway  process  in  simulations  8  and  9,  but  in  these 
cases  the  flame  brush  accelerated  significantly  and  exited  the  domain  before  DDT  could  occur.  The  overall  growth 
rate  of  St  was  lower  than  in  simulation  10  (jea  increases  with  decreasing  [//).  All  simulations  are  well-resolved 
with  resolution  at  least  Ax  =  6^/16  [20].  Convergence  at  this  resolution  was  confirmed  for  the  QSS  in  cases  6 
[20,  21]  and  7  using  resolutions  Ax  =  A /.,o / 8  -  Sl, o/32,  and  convergence  during  the  runaway  was  confirmed  in  case  10 
using  resolutions  Ax  =  6l, 0/8  -  6i„o/ 1 6.  We  also  observed  DDT  in  simulation  4  which,  however,  was  under-resolved 
(Ax  =  <5z,,o/4)  and,  thus,  is  not  shown  here. 

To  determine  the  dependence  of  our  results  on  the  reaction  model,  we  carried  out  a  similar  simulation  in  a  sto¬ 
ichiometric  CFQ-air  mixture.  In  this  case,  d/o  =  0.042  cm  is  close  to  that  in  FF-air,  but  the  flame  speed  is  almost 
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an  order  of  magnitude  lower,  S l,o  =  38  cm/s  [13].  We  also  observed  DDT  in  this  fuel,  but  only  at  higher  turbulent 
intensity  relative  to  S l  with  U\  -  2.24  x  103  cm/s  «  59 Si  and  in  a  larger  system  l  =  0.31  cm  ( L  -  1.328  cm)  (case 
12,  Figs.  1  and  2).  The  overall  evolution,  however,  was  different  from  case  10.  The  time  to  DDT  was  «  and  the 
flame  never  developed  a  QSS.  The  flame  accelerated  significantly  relative  to  fuel,  which  required  a  longer  domain 
(256  x  256  x  8192)  to  observe  DDT,  and,  in  contrast  with  simulation  10,  a  strong  well-defined  global  shock  was 
formed.  This  suggests  that  there  exist  two  distinct  types  of  flame  evolution  in  such  unstable  regimes. 

The  key  aspect  of  the  spontaneous  DDT  mechanism  discussed  here  is  that  it  does  not  place  any  specific  constraints 
on  the  equation  of  state,  reaction  model,  or  the  flame  properties.  Decrease  of  fluid  density  with  increasing  temperature 
in  an  exothermic  process  means  that  at  a  sufficiently  high  but  subsonic  burning  speed,  the  flow  of  products  will  become 
supersonic,  irrespective  of  how  burning  occurs.  This  ensures  that  the  pressure  wave  remains  coupled  to  the  region 
in  which  the  energy  release  occurs  (cf.  location  of  peaks  of  P  and  E  in  Fig.  3b  and  c).  This  is  in  contrast  with  the 
spontaneous  reaction  wave  model  [16]  that  requires  very  specific  hot  spot  properties  in  order  for  the  resulting  reaction 
wave  to  remain  properly  coupled  to  the  pressure  pulse  that  it  generates. 

Figure  1  suggests  that  there  is  both  a  minimal  system  size  and  minimal  relative  turbulent  intensity  at  which  DDT  is 
possible,  and  their  values  are  not  universal  as  they  appear  to  increase  for  reactive  mixtures  with  slower  laminar  flames. 
Applying  the  criterion  given  by  eq.  (1)  to  establish  whether  DDT  can  occur  in  a  given  turbulent  flow  critically  depends 
on  our  ability  to  predict  the  turbulent  flame  speed  for  given  U/  and  /.  This  is  particularly  difficult  in  the  high-speed 
regimes,  and  it  is  in  these  regimes  where  spontaneous  DDT  is  most  likely  to  occur. 
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